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Abstract 

We identify a new and important global (or non- 
binary) constraint. This constraint ensures that the 
values taken by two vectors of variables, when 
viewed as multisets, are ordered. This constraint is 
useful for a number of different applications includ- 
ing breaking symmetry and fuzzy constraint sat- 
isfaction. We propose and implement an efficient 
linear time algorithm for enforcing generalised arc- 
consistency on such a multiset ordering constraint. 
Experimental results on several problem domains 
show considerable promise. 

1 Introduction 

Global (or non-binary) constraints are one of the factors cen- 



tral to the success of constraint programming [Regin, 1994 



Regin,~T996| |Beldiceanu, 2000| . Global constraints specify 



patterns that occur in many problems, and call efficient and 
effective constraint propagation algorithms. In this paper, we 
identify a new and important global constraint. This con- 
straint ensures that the values taken by two vectors of vari- 
ables, when viewed as multisets, are ordered. Such a con- 
straint is useful in a number of domains. For example, in the 
progressive party problem (prob013 in csplib.org), we wish 
to assign a host for each guest and period. We can model this 
with a vector of variables for each period. Each variable is 
assigned the host for a particular guest. This model has un- 
necessary symmetry as the periods can be freely permuted. 
We can break this symmetry by considering the multiset of 
values associated with each vector and ordering these multi- 
sets. The aim of this paper is to study such multiset ordering 
constraints and to develop efficient and effective techniques 
for enforcing them. 

2 Formal Background 

A constraint satisfaction problem (CSP) consists of a set of 
variables, each with a finite domain of values, and a set of 
constraints that specify allowed values for subsets of vari- 
ables. A solution is an assignment of values to the variables 
satisfying the constraints. To find such solutions, we explore 
partial assignments enforcing a local consistency like gener- 
alized arc-consistency (GAC). A constraint is GAC iff, when 



a variable in the constraint is assigned a value, compatible 
values exist for all the other variables in the constraint. GAC 
reduces to arc-consistency (AC) for binary constraints and to 
node-consistency (NC) for unary constraints. Finally, a con- 
straint is bounds consistent (BC) iff, when a variable in the 
constraint is assigned its maximum or minimum value, there 
exist compatible values for all the other variables in the con- 
straint. If a constraint c is NC, BC, AC or GAC then we write 
NC(c), BC(c), AC(c) or GAC(c) respectively. 



We now define the multiset ordering. A multiset is an un- 
ordered list in which repetition is allowed. A multiset M is 
ordered smaller than another N, written M -< m N iff either 
M is empty and N is not, or the largest value in M is smaller 
than the largest value in N, or the largest values are the same 
and, if we eliminate one occurrence of the largest value from 
both M and N, the resulting two multisets are ordered. We 
can weaken the ordering to include multiset equality. That 
is, M ^ m N iff M = N or M -< m N. As in the intro- 
ductory example, we will often view the values taken by a 
vector of variables as a multiset. Given two vectors, x and 
y, we write a multiset ordering constraint as x < m y. This 
constraint is satisfied iff the values for the variables in the 
vectors, when treated as multisets, satisfy the multiset order- 
ing. Similarly, we write a strict multiset ordering constraint as 
x < m V- Throughout the paper, we assume that the variables 
being ordered are disjoint and not repeated. 



We also need the following notation. Input vectors are in- 
dexed from 0. The minimum element in the domain of Xi is 
min(xi), and the maximum is max(xi). The floor function, 
floor(x) assigns all variables in x to their minimum values, 
whilst the ceiling function, ceil(x) assigns all to their maxi- 
mums. The vector x v= d is identical to x except v now has the 
domain {d}. The function occ(x) computes the occurrence 
vector associated with x. We index occurrence vectors in de- 
creasing order from the maximum to the minimum value from 
the domains in x. When comparing two occurrence vectors, 
we assume they start with the occurrence of the same value, 
adding leading zeroes as necessary. Finally, x <i ex y iff x is 
lexicographically less than or equal to y. 



3 Motivating applications 
3.1 Matrix symmetry 

Many constraints programs contain matrices of decision vari- 
ables (so called "matrix models"), and the rows and/or 
columns of these matrices are symmetric and can be per- 
muted flFlener et al, 2002] . Such symmetries are very dif- 
ficult to deal with as there are a super-exponential number 
of permutations of the rows or columns to consider. There 
are several ways to break symmetry in a CSP, such as SBDS 
| |Gent and Smith, 2000| or SBDD [ |Fahle ef a/., 2001] . One 
of the most effective, and the one which we will concentrate 
on as a major application for a multiset ordering constraint, 
is adding extra symmetry-breaking constraints to an initial 
model. Existing techniques for dealing with such symmetries 
typically eliminate only some of the symmetry. Additional 
techniques, like those proposed here, are therefore of consid- 
erable value. 

The progressive party problem mentioned earlier has a 2d 
matrix of decision variables with matrix symmetry. The rows 
of the matrix are the guests, the columns are the periods. Each 
variable gives the host assigned to a given guest in a given pe- 
riod. As periods are indistinguishable, the columns of the ma- 
trix are symmetric. One way to break such column symmetry 
is to lex order the columns | |Frisch et al, 2002") . Similarly, as 
guests can be indistinguishable, (some of) the rows may be 
symmetric and can be lex ordered. Alternatively, we can treat 
each row and/or column as a multiset and break such symme- 
try by multiset ordering the rows and/or columns. 

Unlike lex ordering, multiset ordering the rows of a matrix 
may not eliminate all row symmetry. For example, consider 
the symmetric matrices: 
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Both satisfy the constraint that the first row is multiset less 
than the second. It is therefore a little surprising to dis- 
cover that multiset ordering (which does not break all row 
symmetry) is not dominated by lex ordering (which does) 
but is incomparable. For example, (0,2) <i ex (1,1) but 
(1,1) < m (0,2). 

When we have both row and column symmetry, we can 
multiset order both rows and columns. Like lex ordering both 
rows and columns, this may not eliminate all row and column 
symmetry. Consider the symmetric matrices: 

(: s :) (s : :) 

Both have multiset ordered rows and columns. Unsurpris- 
ingly, multiset ordering rows and columns is incomparable 
to lex ordering rows and columns. Consider the symmetric 
matrices: 
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The first has lex ordered rows and columns, but the columns 
are not multiset ordered. The second has rows and columns 
that are multiset ordered but the columns are not lex ordered. 

An alternative way to deal with row and column symmetry 
is to multiset order in one dimension and apply the symme- 
try breaking method of our choice in the other dimension. 
This is one of the best features of using multiset ordering 



to break symmetry. It is compatible with any other method 
in the other dimension. For instance, we can multiset or- 
der the rows and lex order the columns. Preliminary results 
in iKiziltan and Smith, 20021 suggest that such a combined 
method is very promising. This combined method does not 
eliminate all symmetry (but it is unlikely that any polynomial 
set of constraints does). Consider the symmetric matrices: 

(! 11) (11 I) 

\ 1 2 J \ 2 1 J 

Both have rows that are multiset ordered, and rows and 
columns that are lex ordered. Multiset ordering the rows and 
lex ordering the columns is again incomparable to lex order- 
ing rows and columns. Consider the symmetric matrices: 
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The first matrix has rows that are multiset ordered and 
columns that are lex ordered. However, its rows are not lex 
ordered. The second matrix has rows and columns that are 
lex ordered but does not have rows that are multiset ordered. 
Whilst the two orderings are theoretically incomparable, our 
experimental results (see later) show that multiset ordering 
the rows and lex ordering the columns is often the most ef- 
fective symmetry breaking constraint currently known. 

3.2 Fuzzy constraints 

Another application for multiset ordering is to fuzzy CSPs. A 
fuzzy constraint associates a degree of satisfaction to an as- 
signment tuple for the variables it constrains. To combine de- 
grees of satisfaction, we can use a combination operator like 
the minimum function. Unfortunately, the minimum function 
may cause a drowning effect when one poorly satisfied con- 
straint 'drowns' many highly satisfied constraints. One solu- 
tion is to collect a vector of degrees of satisfaction, sort these 
values in ascending order and compare them lexicographi- 
cally. This leximin combination operator identifies the as- 
signment that violates the fewest constraints |Fargier, 19941. 
This induces an ordering identical to the multiset ordering ex- 
cept that the lower elements of the satisfaction scale are the 
more significant. It is simple to modify a multiset ordering 
constraint to consider the values in a reverse order. To solve 
such leximin fuzzy CSPs, we can then use branch and bound, 
adding an ordering constraint when we find a solution to en- 
sure that future solutions are greater in the leximin ordering. 

4 GAC algorithm for multiset ordering 

The last section motivated why we want multiset ordering 
constraints. We need, however, to be able to propagate such 
constraints efficiently. We therefore developed an efficient 
GAC algorithm for such constraints. 

4.1 Background 

The algorithm exploits two theoretical results. The first re- 
duces the problem to testing support for upper bounds of x 
and lower bounds of y on suitable ground vectors. The sec- 
ond reduces these tests to lex ordering suitable occurrence 



vectors. Identical results hold for the strict multiset ordering 
constraint but for reasons of space we omit them here. 
Lemma 1 Given two disjoint and non-repeating vectors of 
variables, x and y, with non-empty domains, GAC{x < m y) 
iff Vxi e x, y-j £ y: 

floor(x x . ) < m ceil(y) 

floor(x) < m ceil(y y . 

Proof: (=K) As the constraint is GAC, all values have sup- 
port. In particular, Xi — max(xi) has support. The best 
support comes if all the other variables in x take their mini- 
mums, and all the variables in y take their maximums. Hence, 
floor( K x Xi — raax ^ Xi 'j ) < m ceil(y). Similarly, for y t . 

(4=) The first constraint ensures that max(xi) is supported. 
The values which support max(xi) also support all values 
smaller. Hence, all the values in the domain of Xi are sup- 
ported. By an analogous argument, all the values in the do- 
main of yi are supported. Hence the constraint is GAC. QED. 

The next lemma reduces these test for support to lex order- 
ing suitable occurrence vectors. 

Lemma 2 Given two multisets of values, M and N, M < rn 

N iffocc( M) < iex occ(N). 

Proof: See | |Kiziltan and Walsh, 20021 . 

4.2 A worked example 

Based on these lemmas, we have designed an efficient algo- 
rithm for enforcing GAC on a multiset ordering constraint. 
The algorithm goes through the Xi and yj checking for sup- 
port in the appropriate occurrence vectors. Incremental com- 
putation of the lex ordering test avoids repetition of work. 
Consider the multiset ordering constraint x < m y where: 

3 = <{5},{4,5},{3,4,5},{2,4},{1},{1}) 
y = <{4,5},{4},{1,2,3,4},{2,3},{1},{0}> 

We construct occurrence vectors for floor(x) and ceil(y), in- 
dexed from 5 to 0: 

occ(floor(x)) = (1,1,1,1,2,0) 
occ(ceil(y)) = (1,2,1,0,1,1) 

Recall that index i in occ(ceil(y)) denotes the number of 
occurrences of the value i in ceil(y). For example, index 4 is 
2 as the value 4 occurs twice. 

We first check if occ(floor(x)) >i ex occ(ceil(y)). If so, we 
can fail immediately because no value for any variable can 
have support. Here, occ(floor(x)) <i ex occ(ceil(y)). In fact, 
we record (in a pointer, a) that the two occurrence vectors 
are lex ordered by index 4 of occ(floor(x)), which is strictly 
smaller than index 4 of occ(ceil(y)). This means that we will 
fail to find support in the yj if any of the xi is assigned a new 
value greater than 4. We now go through the Xi checking for 
support for their maximum values, and then the yi checking 
for support for their minimum values. 

Consider xq. As it has a singleton domain, and 
occ(floor(x)) <i ex occ(ceil(y)), its only value must 
have support so we skip it. Now consider x%. Do 
its values have support? Changing occ(floor(x)) to 
occ((min(xo) , max(xi) 1 . . . , min{x n -i))) decreases the 
number of occurrences of mm(xi) = 4 by 1, and increases 
the number of occurrences of max(xi) = 5 by 1. As 



min(xi) > a — 4, this upsets the lex ordering of the two oc- 
currence vectors. We therefore prune all values in the domain 
of X\ larger than a. This leaves a single supported value, 4. 

Now consider x-i. Changing occ(floor(x)) to 
occ({min(xo) , . . . , max{x2), ■ ■ ■ , min{x n -i))) decreases 
the number of occurrences of min{x2) = 3 by 1, and 
increases the number of occurrences of max(x2) = 5 
by 1. As with xi, any value of x-i larger than a = 4 
upsets the lex ordering. We therefore prune 5 from the 
domain of X2- Now consider X3. Changing occ(floor(x)) to 
occ{{min(xo), . . . , max{x^), . . . , min{x n -i))) decreases 
the number of occurrences of min{xz) — 2 by 1, and 
increases the number of occurrences of max{x^) = 4 by 1. 
The occurrence vectors beneath a would now be lex ordered 
the wrong way. We therefore also prune the value a = 4, 
leaving a single supported value 2 in the domain of 23. As 
X4 and £5 have singleton domains, their values have support. 

Similarly, we check the minimums of the yj for support. 
However, rather than prune values above (and in some cases 
equal to) a, there is now a dual pointer j3 and we prune values 
in the domains of yj up to (and in some cases equal to) (3. The 
pointer (3 is the largest index such that the occurrence vectors 
beneath it are lex ordered the wrong way. Any value less than 
(3 cannot hope to change the lex ordering as the value at j3 
will still order the vectors the wrong way. Such values can 
therefore be pruned. Once we have considered each of the yj, 
we have the following generalized arc-consistent vectors: 

x = ({5},{4},{3,4},{2},{1},{1}) 
V = <{5},{4},{3,4},{2,3},{1}.{0}) 

4.3 Algorithm details 

The algorithm uses two pointers a and (3, and two flags 7 
and 8 to avoid traversing the occurrence vectors each time we 
look for support. The pointer a is set to to the most significant 
index above which all occurrences are pair-wise equal and at 
a they are strictly ordered. If the vectors are equal then a is 
set to —00. The pointer (3 is set to the most significant index 
below a such that the occurrence vectors are lex ordered the 
wrong way. If no such index exists, we set (3 to —00. The 
flag 7 is set to true if all the indices between a and (3 are 
pair-wise equal and the flag i5 is set to true if the sub-vectors 
below (3 are lex ordered the wrong way. For example, given 
the occurrence vectors in section |4T2l a is set to 4, (3 to 2, and 
the flags 7 and 8 are set to true. 

We summarise the major steps the algorithm performs: 

A. Build occ(floor(x)) and occ(ceil(y)) 

B. Set a, (3, 7, 8 according to their definitions 

C. For each xc If its maximum disturbs the lex ordering 

on the occurrence vectors, tighten its upper-bound to a 
when the occurrence vectors are lex ordered below a, 
otherwise to a — 1. 

D. For each yf. If its minimum disturbs the lex ordering on 

the occurrence vectors, then tighten its lower-bound to 
(3 when the occurrence vectors are lex ordered below a, 
otherwise to f3 + 1. 

When we prune a value, we do not need to check re- 
cursively for previous support. Pruning changes neither the 
lower bounds of x nor the upper bounds of y. These values 



continue to provide support. The exception is when a domain 
is a singleton, and pruning causes a domain wipe-out. 

We now give pseudo-code for an algorithm that maintains 
GAC on a multiset ordering constraint between vectors x and 
y which are of length n and m respectively. As the algorithm 
reasons about occurrences vectors, the original vectors need 
not be identical in length (though they often are). 

The algorithm is called whenever lower bounds of Xi or 
upper bounds of yj change. Lines Al to A3 build the oc- 
currence vectors ox and 6y. Line Bl calls the procedure to 
set the pointers a and j3, and the flags 7 and 5. Lines Cl- 
13 check support for the maximums of the Xi's while lines 
Dl-14 check support for the minimums of the y^'s. 

Procedure GACMSOO 

Al. I :— min({min(xi)\i G [0, n)} U {min(yj)\j G [0, m)}) 

A2. u :— max({max(xi)\i G [0, n)} U {max(yj)\j G [0, m)}) 

A3, ox — occifioor(x)), 6y — occ(ceil(y)) 

Bl. SetPointers-and-Flags(i, ti) 

CI. FOR i = TO 71 — 1 % check support for x's 

C2. IF (min(xi) ^ max(xi)) 

C3. a :— min(xi) b :— max(xi) 

C4. IF (a > a) NC(xi = a) 

C5. IF (a > a > /3 & 6 > a) NC(x, < a) 



C6. 


IF (a = & b > a) 


C7. 


IF (ox a + 1 — oy a ) 


C8. 


IF (7 & oxf, - 1 = oyp & S) NC(xi < a) 


C9. 


ELSE NC(x, < a) 


CIO. 


IF (a < & b > a) 


Cll. 


IF(ox a + 1 = oy a ) 


C12. 


IF (7) NC(xt < a) 


C13. 


ELSE NC(xi < a) 


Dl. 


FOR j — TO m — 1 % check support for y's 


D2. 


IF (min(yj) ^ max(yj)) 


D3. 


a :— min(yj) b :— max(yj) 


D4. 


IF (6 > a)NC( Vj = b) 


D5. 


IF (fe = a & oz Q + 1 = ojf Q ) 


D6. 


IF(a < /3 & 7) 


D7. 


IF (ox ,3 = oy,, + 1) 


D8. 


IF(5)JVC( % > /3) 


D9. 


ELSE NC(y 3 > 0) 


D10. 


ELSE JVC( % > 0) 



Procedure SetPointers-and-Flags(i, u) 

Bl. 7 :=false, 5 :=false, a :— u 

B2. WHILE (a>!A oa; Q = oy a ) a := cv - 1 

B3. IF (a > / A oi a > oy a ) FAIL 

B4. IF (a = / - 1) a := -oo, /3 := -co 

B5. ELSE 

B6. := a — 1, 7 -.—true 
B7. WHILE (ft > i A 01,3 < o W ) 
B8. IF {ox p < oyp) 7 : = false 
B9. /3 := /3 - 1 

BIO. IF (/3 = / - 1) /3 := -oo, 7 : = false 
Bll. IF (/3 # -00 A0> I) 
B12. i := /3 - 1 

B13. WHILE (i > ! A 01, = o^) i := i - 1 
B14. IF (i > I A oo;; > oi/i) 5 := Jrae 

For each x i7 we only check for support if Xi is not singleton 
(line C2). There are six cases where we prune the domain of 
XC (1) at line C4, as min(xi) > a, any value in the domain 
of Xi greater than min(xi) lacks support because it will dis- 
turb the lex ordering; (2) at line C5, as a > min(xi) > /?, 
and max(xi) > a, all the values in the domain of Xi greater 



than a disturb the lex ordering, and lack support; (3) at line 
C6, C7, and C8, all values greater than or equal to a lack 
support. If we assign a to x i7 then the vectors will be equal 
at index a and (3, the values between a and [3 are also all 
pair-wise equal (since 7 is true), and the vectors below (3 are 
ordered the wrong way (since 5 is true). Thus, the value a 
also lacks support and is pruned; (4) at line C9, all values 
greater than a lack support; (5) at line C10, Cll, and C12, 
all values greater than or equal to a lack support. If we assign 
a to the vectors will be lex ordered the wrong way as the 
difference between the number of occurrences at a is exactly 
one, and 7 is true. Thus, the value a also lacks support and is 
pruned; (6) at line C13 all values greater than a lack support. 

For each yj, we only check for support if yj is not singleton 
(line D2). There are four cases where we prune the domain 
of yf (1) at line D4, as max(yj) > a, any value smaller than 
max(yj) disturbs the lex ordering and lacks support; (2) at 
lines D5 to D8, the situation is dual to the third case for the 
Xi, and any value less than or equal to f3 lacks support; (3) at 
line D9, any value less than (3 lacks support; (4) at line D10, 
again, any value less than or equal to (3 lacks support. 

4.4 Theoretical properties 

In a longer technical report, we prove the following results 
about the algorithm's complexity and correctness. 

Theorem 1 GACMSO runs in time 0(n + m + d) where d 
is u — I. 

If d n,m (and for multisets, we expect this as the num- 
ber of values is typically less than the cardinalities to permit 
repetition), then the algorithm is 0(n + m). 

Theorem 2 For disjoint and non-repeating vectors, 
GACMSO either establishes failure if x < m y is not 
satisfiable, or prunes values from x and y to ensure 
GAC{x < m y). 

The algorithm can easily be modified to support strict mul- 
tiset ordering. The only differences are that we fail if ox = 6y 
and that (3 is set to I — 1 not —00. The algorithm then en- 
forces a strict inequality on the occurrence sub-vectors above 
(3. Another variant of the algorithm is when d » n,m. In 
such a situation, it could be costly to construct the occurrence 
vectors. We can instead sort the minimums of the Xi and 
the maximums of the yj, and compute a, (3, 7 and 5 as if 
we had the occurrences by scanning these sorted lists. This 
information is all we need to compute support for each Xi 
and yj in turn. The complexity of this modified algorithm is 
0(n log n + m log m) as the cost of sorting dominates. 

If we have multiple rows of a matrix that are multiset or- 
dered, we can decompose this into multiset ordering con- 
straint on all pairs of rows, or (further still) onto ordering 
constraints just on neighbouring pairs of rows. The following 
result shows that such decompositions hinder constraint prop- 
agation. Nevertheless, it will usually be most cost effective to 
post just the 0{n) ordering constraints between neighbouring 
pairs rather than the 0(n 2 ) constraints between all pairs. 

Theorem 3 GAC(ii < j . x t < m Xj) is strictly stronger 
than GAC(x*i < m Xj ) for all i < j, and this itself is strictly 
stronger than GAC(xi < m x i+ i) for all i. 



5 Alternative approaches 

5.1 Arithmetic constraint 

Barbara Smith (personal communication) has proposed en- 
forcing x < m y on vectors of length n via the arithmetic 
constraint n x " + . . . n Xn - 1 < n v " + . . . rt^"- 1 . This is similar 
to the transformation of a leximin fuzzy CSP into an equiva- 
lent MAX CSP |Schiexef aT7l 995 |. BC on such a constraint 
is equivalent to GAC on the original multiset ordering con- 
straint. However, such an arithmetic constraint is only feasi- 
ble for small n. Further, most existing solvers will not enforce 
BC on such an arithmetic constraint, but will delay it until all 
but one of the variables are instantiated. 

5.2 Decomposition 

Multiset ordering is equivalent to the lex ordering the as- 
sociated occurrence vectors. As we have efficient algo- 
rithms for constructing occurrence vectors (via the global 
cardinality constraint jRegin, 1996| ) and for lex ordering 
| |Frisch et ai, 2002) , this might be an alternative approach. 
However, as the following theorem shows, such a decomposi- 
tion hinders constraint propagation. Also, the two global car- 
dinality constraints in such a decomposition are more expen- 
sive to enforce than the algorithm presented here. We write 
gcc(x, ox) for the global cardinality constraint that channels 
from a vector of variables x to the associated occurrence vec- 
tor ox. 

Theorem 4 GAC(x < rn y) is strictly stronger than simulta- 
neously enforcing GAC(gcc(x, ox)), GAC(gcc(y, dy)), and 
GAC(6x < lex 6y). 

Proof: Clearly it is as strong. To show strict- 
ness, consider x = ({1, 2}, {1, 2}, {2}, {2}) and y = 
({1, 2}, {1, 2}, {0, 1, 2}, {0, 1}}. The multiset ordering con- 
straint is not GAC since in yi has no support but the decom- 
position is unable to determine this. QED. 

Another approach is to use the sorted constraint in the 
Eclipse solver. This ensures that the values taken by one vec- 
tor of variables are identical but in sorted order to the values 
taken by a second vector of variables. To post a multiset or- 
dering constraint on two vectors, we can channel each into a 
sorted vector and lex order these. The above example demon- 
strates that such a decomposition again hinders propagation. 
The sorting constraint is also more expensive to enforce. 

6 Experimental results 

We designed some experiments to test three goals. First, is 
multiset ordering an effective method for dealing with row 
and/or column symmetry? Second, how does multiset order- 
ing compare to lex ordering? Which one breaks more sym- 
metry? Is a combined method, which multiset orders one di- 
mension and lex orders the other one of the matrix, superior? 
Third, does our GAC algorithm do more inference in prac- 
tice than the decomposition? Similarly, is the algorithm more 
efficient in practice than its decomposition? 

The symmetry breaking constraints we used are strict lex 
ordering on the columns (<i ex C), on the rows (<i ex R); mul- 
tiset ordering on the rows (< m R), (strict) multiset ordering 
on the columns (< m C and < m C); and combinations of these 
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Table 1: 5-13-29 progressive party problem using row -by- 
row labelling. 

constraints. Such constraints are posted between adjacent 
rows/columns. The results of the experiments are shown in 
tables where a "-" means no result is obtained in 1 hour (3600 
sees). The experiments are done using ILOG Solver 5.2 on a 
1000MHz pentimum III with 256 Mb RAM. 

6.1 Progressive Party Problem 

There are a set of host boats, each with a capacity, and a set of 
guest boats, each with a crew size. We wish to assign a host 
for each guest and period, such that a guest crew never visits 
the same host twice, no two guest crews meet more than once, 
and the spare capacity of each host boat, after accommodating 
its own crew, is not exceeded (prob013 in csplib.org). 

A matrix model of this problem | |Smith et ai, 1995] is a 
2-d matrix of guests x periods where each variable is as- 
signed a host representing that a host is accommodating a 
particular guest in a given time period. The rows are the 
guests, the columns are the periods. This model has col- 
umn and partial row symmetry: any two periods, and any 
two guests with the same crew size are indistinguishable. We 
consider the 13-hosts and 29 guests problem with 5 and 6 
time periods, referred as 5-13-29 and 6-13-29. These prob- 
lems have p!l4!2!4!5!7! row and column symmetries where 
p is the number of time periods. The actual data can be 
found in csplib.org. Due to the problem constraints, no pair 
of rows/columns can be equal, hence we can safely pose strict 
lex ordering. However, any two distinct rows/columns might 
be equal when viewed as multisets. 

As in [S mith et ah, 1995| , the guest boats are ordered in 
descending order of their size. We order the host boats in 
descending order of spare capacity to choose a value in a 
succeed-first manner. Results obtained by row-by-row, and 
column-by-column labelling strategies are given in Tables Q] 
and |2 With row -by-row labelling, we cannot solve 6-13-29 
with or without symmetry breaking. For the 5- 1 3-29 problem, 
<i ex R breaks a lot more row symmetry than < m R. However, 
the reverse is true for the columns. Here, <i ex C does not 
break any symmetry but < m C does. Multiset ordering one 
dimension of a matrix therefore does not necessarily break 
less symmetry than lex ordering the same dimension. Such 
phenomena occur through interactions with the search strat- 
egy: a search strategy might already lex order, so multiset 
ordering constraints break additional symmetry. The smallest 
search tree and also the least solving time is obtained by < m C 
<; ex R. This supports our conjecture that lex ordering one di- 
mension combined with multiset ordering the other can break 
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Table 2: 6-13-29 progressive party problem using column- 
by-column labelling. 
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Table 3: Sports scheduling problem. 



more symmetry than lex/multiset ordering both dimensions. 

With column-by-column labelling, we are able to solve the 
6-13-29 problem. Neither of </ ea; R, <i ex C, <i ex RC break 
any symmetry. The smallest search tree is obtained by < m C. 
This supports our conjecture that multiset ordering one di- 
mension can break more symmetry than lex ordering the same 
or both dimensions. If the search strategy already orders both 
dimensions lexicographically, imposing a constraint like mul- 
tiset ordering in one dimension breaks additional symmetry. 

6.2 Sports Scheduling with Odd Teams 

This is a modified version of prob026 in csplib.org. We have 
n teams (n is odd), playing over n weeks. Each week is 
divided into (n — l)/2 periods, and each period is divided 
into 2 slots, home and away. We wish to find a schedule 
so that every team plays at most once a week, every team 
plays twice in the same period over the tournament and every 
team plays every other team. We slightly modify the model 
in I Van Hentenryck et ai, 19991, where teams is a 3-d matrix 
of periods x weeks x slots. Each element of teams is the 
team playing in a given period, week and slot. We treat this 
matrix as 2-d where the rows are the periods and columns are 
the weeks, and each entry is a list of variables giving the slots. 

As the periods and the weeks are indistinguishable, this 
problem has n\(n — 1/2)! row and column symmetries. We 
pose strict ordering constraints on the rows and columns of 
teams as the periods and weeks cannot be equal. Due to the 
constraints on the periods, posing multiset ordering on the 
rows is not effective. 

Results obtained by column-by-column labelling of the 
teams are given in Table [3] For one column, we first label 
the first slots; for the other, we first label the second slots. 
With this strategy, </ ea; R does not break any symmetry, so 
we omit it in the table. Posing multiset ordering by our algo- 
rithm is much more effective and efficient than by gec and lex 
ordering constraints. This holds for many other search strate- 
gies. In Table[3] we note that < m C gives a smaller search tree 
than <i ex C. However, for other search strategies the reverse 
is true. This supports the theoretical result that lex ordering 
and multiset ordering are incomparable. 

7 Conclusions 

We have identified a new and important global (non-binary) 
constraint. This constraint ensures that the values taken by 
two vectors of variables, when viewed as multisets, are or- 
dered. We have developed an efficient linear time algorithm 



for enforcing generalised arc-consistency on such a multiset 
ordering constraint. We have proposed a number of appli- 
cations for this new constraint including breaking symmetry 
in matrix models, and fuzzy constraint satisfaction. We have 
shown that alternative methods for posting a multiset ordering 
constraint like an arithmetic constraint or decomposition are 
inferior. Finally, we have implemented this generalized arc- 
consistency algorithm in ILOG Solver. Experimental results 
on a number of problem domains show considerable promise. 
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